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6.4 Geothermal convection in porous media 

Motivations: Geothermal Energy, Pipe hne in permafrost, thermal stress around buried 
nuclear waste. 

S. Chandrasekhar: Hydrodynamic and Hydromagnetic Stabililty, Oxford, 1961. 
D. A. Nield and A. Bejan, Convection in porous media , Springer- Verlag, 1992. 
Let us begin with the effective equations over macroscales. 
Recall the definitions of seepage velocity, 

/// ^dV (6.4.1) 

and of porosity: 



n = ^ (6.4.2) 
Mass conservation for an incompressible fluid: 

V-u = (6.4.3) 

Momentum conservation (Darcy 's law) 

= -Vp-^u + p/g (6.4.4) 

Inview of the small size of the solid grains heat diffuses by conduction almost instantly in 
time scales of geothermal interest. We therefore assume the temperature in both phases to 
be equal, 

Ts^Tf^T (6.4.5) 

The energy conservation law in each phase is : 
Solid 

cfT 

(1 - n){pC)s^ = (1 - n)KsV^T (6.4.6) 

fluid: 

dT 

<pC)f-^ + {pC)fn ■ VT = nKfV^T (6.4.7) 
Adding the two we get the energy equation for the mixture 

(pC7)^— + {pC) fU-VT = K^V^T (6.4.8) 
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where 

{pC)m = (1 - n){pC)s + n{pC) f, Km^{l- n)Ks + nKf (6.4.9) 

Defining 



ipC)f 

and 



(6.4.10) 



«„=(^ (6.4.11) 

as the mixture diffusivity, we can rewrite (6.4.8) as 

dT 

(7— + u • VT = K^V^T (6.4.12) 
at 

Note: Strictly speaking, the above equation ought to be derived from the microscale 
considerations, i.e., by homogenization analysis. 
Equation of state: 

Pf^Po[l-f3{T-T„)] (6.4.13) 
In summary the governing eqautions are (6.4.3), (6.4.4), (6.4.12) and 6.4.13). 



